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ABSTRACT 



In this paper we provide a possible link between the structure of the bipolar nebula 
surrounding SN1987A and the properties of its progenitor star. A Wind Blown Bubble 
(WBB) scenario is employed, in which a fast, tenuous wind from a Blue Supergiant 
expands into a slow, dense wind, expelled during an earlier Red Supergiant phase. 
The bipolar shape develops due to a pole-to-equator density contrast in the slow wind 
(i.e., the slow wind forms a slow torus). We use the Wind Compressed Disk (WCD) 
model of Bjorkman & Cassinelli (1992) to determine the shape of the slow torus. 
In the WCD scenario, the shape of the torus is determined by the rotation of the 
progenitor star. We then use a self-similar semi-analytical method for wind blown 
bubble evolution to determine the shape of the resulting bipolar nebula. 

We find that the union of wind-compressed-disk and bipolar-wind-blown-bubble 
models allows us to recover the salient properties of SN1987A's circumstellar nebula. 
In particular, the size, speed and density of SN1987A's inner ring are easily reproduced 
in our calculations. An exploration of parameter space shows that the red supergiant 
progenitor must have been rotating at ^ 0.3 of its breakup speed. We conclude that 
the progenitor was most likely spun up by a merger with a binary companion. Using a 
simple model for the binary merger we find that the companion is likely to have had a 
mass ;> O.5M . 
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1. Introduction 

The rings surrounding SN1987A (Burrows et al. 1995) are one example of the ubiquitous 
phenomenon known as a bipolar outflow. These hypersonic, figure-8 shaped nebulae occur in 
almost all forms of evolved stars, regardless of their mass {e.g., Planetary Nebulae: Manchado 
et al. 1996; Schwartz, Corradi & Melnick 1992; Balick 1987, Luminous Blue Variables: Nota et 
al. 1995; and see Livio 1997 for a review). 

Considerable progress has been made in understanding bipolar outflows resulting from 
interacting stellar winds. Analytical and numerical studies have recovered many of the observed 
features of these outflows through the Generalized Wind Blown Bubble (GWBB) paradigm (Kwok 
& West 1984; see Frank 1998 for a review). In this scenario, a fast wind from a central source 
expands into a highly aspherical (i.e., toroidal) environment. The interaction of the wind and its 
environment produces an expanding bubble bounded by strong shocks. The bubble's velocity is 
highest in the direction of lowest density. As a result, the density gradient in the environment 
(the slow torus), establishes a preferred axis for the bipolar lobes. The GWBB paradigm has 
been successfully applied to all forms of evolved-star bipolar outflows. Models which include 
the relevant hydrodynamics and microphysics have recovered the global morphology, kinematics 
and ionization patterns in many planetary nebulae (PNe; Frank & Mellema 1994, Mellema 1996, 
Dwarkadas, Chevalier & Blondin 1995). Clear correspondences also exist between GWBB models 
and the shapes and kinematics of Wolf-Rayet (WR) nebulae (Garcia-Segura & MacLow 1993), 
LBVs like rj Carina (Frank, Balick & Davidson 1995) and symbiotic stars like R Aquarii (Henney 
k Dyson 1992). 

The success of the GWBB paradigm opens the possibility of using the properties of a nebula 
to infer the history of the central star in terms of its mass loss. We note that as of yet no studies 
have attempted to make an explicit and quantitaive link between the properties of the nebula 
and history of the star. The intense scrutiny applied to SN87A makes it a unique laboratory for 
studying the connection between bipolar outflows and stellar evolution. Such is the goal of this 
paper. The recent discovery of a bipolar outflow surrounding Sher 25 (Brander et al. 1997), a 
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star similar to the progenitor of SN87A, also raises the possibility that SN87A is not an isolated 
object, but defines a new class of bipolar outflows. 

When the central ring of SN87A was discovered, it was quickly interpreted as the waist of a 
bipolar wind-blown bubble (Luo & McCray 1991, Wang & Mazzuli 1992; a structure anticipated 
by Soker h Livio 1989). In these models the ambient medium was taken to be a toroidal wind 
deposited by SN87A's progenitor in its Red Supergiant (RSG) stage. The bubble was subsequently 
inflated by a fast wind from the star's penultimate incarnation as a Blue Supergiant (BSG). The 
first numerical simulations of this process were carried out by Blondin & Lundqvist (1993: BL93) 
who used an ad hoc function to determine the asphericity of the slow RSG wind. BL93 were 
quite successful in demonstrating the feasibility of the GWBB paradigm for SN87A. However, the 
observed low expansion speed of the central ring (v m ~ 8 km s _1 , Meaburn et al. 1995) presented 
a problem. To create a model with the correct kinematics, BL93 were forced to adopt very low 
values of both the RSG and BSG wind velocities. In particular, their value of the RSG wind 
speed, v a ~ 5 km s _1 , seemed particularly anomalous. Canonical values for Red Giant winds are 
comparable to 20 km/s (Habing 1996). BL93 were also forced to take an equator-to-pole density 
contrast of p = p e /pp = 20, which, at the time, seemed large. The size of p was one reason cited 
by McCray & Lin 1992 in their arguments that the ambient density distribution represented a 
remnant protostellar disk rather than a stellar wind (although more recently the protostellar disk 
idea has been abandoned, McCray, private communication). Martin and Arnett (1994, MA95) 
carried out GWBB simulations similar to Blondin & Lundqvist. Using a different ad hoc function 
to control the density asymmetry of the RSG wind, MA95 confirmed Blondin & Lundqvist's 
results, including the need for low RSG wind speeds, v a < 10 km/s. We note that both BL93 and 
MA95 assumed that the the RSG wind speed was isotropic. 

The discovery of additional upper and lower rings of SN87A both confirmed and confused 
the image of SN87A's nebula as a bipolar outflow. While a number of different models for 
the the formation of these outer rings have been proposed (e.g., Burrows et al. , 1995; MA95; 
Podsiadlowski, Fabian & Stevens, 1991; Lloyd, O'Brian &: Kahn, 1995; Meyer, 1997), there does 
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not yet exist a generally accepted model. 

In spite of the difficulties associated with the upper and lower rings, the consensus appears 
to be that the GWBB model has proven effective for SN87A (Crotts 1998). Its success has 
allowed researchers to pose basic questions about the progenitor. In particular, the need for an 
aspherical RSG wind has led a number of authors to posit the existence of a binary companion 
for the progenitor of SN87A. Soker (1998) and Livio (1998) have argued that binary companions 
are required to produce slow toroids in PNe; similar arguments hold for SN87A (Livio 1997). 
Using considerations from stellar evolution theory, along with the need for a high equator to pole 
contrast, Podsiadlowski (1992) also argued that SN87A was a binary system. He calculated a likely 
companion mass of 3 to 6 M with the progenitor /primary of 16 M . Podsiadlowski concluded 
that the two stars merged before the supernova, which is consistent with the lack of observational 
evidence for a companion in the post-supernova epoch (e.g., Crotts, Kunkle & Heathcote, 1995). 

In this paper, we take up the challenge posed by previous models: whereas previous 
investigations have used an ad hoc description of the aspherical RSG, we will link the shape of this 
slow torus to the properties of the star itself. Currently, the best generic scenarios for developing 
asphericity in a slow wind are Common Envelope interactions in binary stars (Livio & Soker 1988) 
and the Wind Compressed Disk (WCD) model of Bjorkman & Cassinelli (1993). In the Common 
Envelope model, a slow, dense wind which has higher mass loss in the orbital plane of the binary 
system is ejected by the primary. While this process has been shown to be effective (Terman et 
al. , 1995; Rasio & Livio, 1996), calculations of the wind shape require numerical models which 
must span a large range in both length and time scales. Thus, there is no simple means for linking 
initial states of the binary with the shape of the ejected slow torus. The WCD model, which was 
developed for excretion disks surrounding B[e] stars, relies on the equatorial focusing of wind 
streamlines from rapidly rotating stars. The advantage of the WCD mechanism is its relatively 
simple formalism, which relates the properties of a single star (mass, temperature, mass loss rate, 
rotation rate) with the properties of the slow torus. The WCD mechanism has been adapted to 
Red Giant stars (Ignace, Bjorkman & Cassinelli, 1996) and has already been used to produce slow 
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tori in LBV and PNe bipolar outflow simulations (Garcia-Segura et al. , 1997a; Garcia-Segura et 
al. , 1997b). 

The goal of the present paper is to utilize the WCD and GWBB formalisms to establish the 
efficacy of the the combined model for SN87A. We then work backwards to bracket the properties 
of the progenitor based on the comparison of models with observations. In what follows, we 
describe self-similar models of the nebula surrounding SN87A, which employ a WCD model to 
determine the properties of the RSG wind. Our self-similar model is based on the work of Giuliani 
(1982) and Dwarkadas, Chevalier & Blondin (1995; hereafter DCB95). We use observations of 
the rings surrounding SN87A to constrain the parameters of the system, and to explore the 
implications for the rotation rate of the progenitor. In §2 we describe the WBB and WCD models. 
In §3 the results of the models and a comparison with observations are presented, as well as 
a "fiducial," or "best-fit" model for the nebula. In §4 the implications of these results for the 
progenitor are explored. Lastly, in §5 we discuss our results. 



2. Theory 

When a stellar wind "turns on," it expands ballistically until enough ambient material is swept 
up for significant momentum to be exchanged between the wind and the ambient medium (Koo 
& McKee 1992). A triplet of hydrodynamic discontinuities then form, defining an "interaction 
region" bounded internally (externally) by undisturbed wind (ambient) gas. The outer boundary 
is an outward-facing shock. It accelerates, compresses, and heats the ambient material as it 
propagates. We refer to this as the ambient shock, and denote its position as R. The inner 
boundary is defined by an inward-facing shock which decelerates, compresses and heats the stellar 
wind. We refer to this feature as the wind shock. Its position is R ws . A contact discontinuity 
(CD), R c d, separates the shocked wind and shocked ambient material. In the 1 dimensional (1-D) 
bubble these discontinuities form a sequence in radius: Rws < Red < R- 

The compressed gas behind either or both shocks emits strongly in optical, UV and IR 
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wavelengths producing a bright shell which defines the observable "bubble." In SN87A, compressed 
material at the equator defines the central ring. The dynamics of the bubble and its emission 
characteristics are defined by the efficiency of post-shock shock cooling. Behind each shock we 
can define a cooling time-scale t c = E t /E t , where E t is the thermal energy density of the gas. 
Radiative cooling can be expressed in terms of a cooling function: E t = C(T) = n 2 A(T), where n 
is the number density of the gas, T its temperature, and the function A(T) is a sum over many 
radiative processes emitting at a variety of wavelengths. The bubble has cooling and dynamical 
time scales defined as 

3kT 

tc ~ 2nT(T) (1) 
U = £, (2) 

where V is the speed of the ambient shock, k is Boltzmann's constant, and we have assumed (as 
we shall throughout this paper) that 7 = 5/3. Comparison of t c and td separates WBBs into two 
classes: Radiative (referred to as momentum conserving) and Adiabatic (or energy conserving). 

In what follows we assume that the BSG wind speed is large enough that the shocked wind 
is extremely hot, and cools inefficiently i.e., t c 3> t^. In this case the wind shock is adiabatic 
and hot shocked wind material fills the expanding bipolar bubble, causing expansion through 
thermal pressure, so that a hot bubble forms. We also assume that the shocked ambient material 
cools efficiently, so that R c d ~ R. DCB95 showed that for thermal-pressure-supported, or 
energy conserving bubbles, the expansion is self-similar for a significant portion of its evolution. 
The resulting, simplified equations may then be solved without resorting to time-dependent 
hydrodynamical calculations. 

The parameters of this wind-blown bubble (WBB) model include the description of the 
BSG and RSG winds. Specifically, the wind speeds and mass loss rates of the two winds must 
be supplied, as well as the dependence of the wind speeds on the polar angle 9, given by the 
asymmetry functions f(9) for the RSG wind density, and g(9) for the RSG wind speed, defined by 

g(9) = lim v r /v a , (3) 
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f(9) = Km 4irr 2 v aP /M a . (4) 

In previous models of SN87A, these asymmetry functions have been determined ad hoc (e.g., Luo 
& McCray, 1991; BL93). Here we assume that the balance of wind driving and centrifugal forces 
focuses the RSG wind toward the equatorial plane, creating a wind-compressed disk (WCD: 
Bjorkman and Cassinelli, 1993). A solution of the WCD equations yields f(9) and g(9). We note 
again that all previous studies assumed an isotropic RSG wind velocity, so that g(9) is constant, 
whereas the WCD model naturally proscribes a wind asphericity. In the following sections we 
describe the models used for the wind blown bubble and the environment. 



2.1. Wind-Blown Bubble Model 

The general equations describing the evolution of a thin shell have been derived by Giuliani 
(1982). In the DCB95 model, these equations are applied to a shell of swept-up, shocked 
Asymptotic Giant Branch star wind material which forms the bright optical regions of Planetary 
Nebulae. Here we apply the model to the swept-up RSG gas. It is assumed that the width of the 
shell is much smaller than its radius, and that magnetic fields may be neglected. The velocity and 
density are averaged over the width of the shell. We briefly sketch the derivation of the equations 
here; the complete derivation may be found in DCB95. 

The continuity equation for the shell is given by 

do . . . . <91n^4 1 d / . „ , ,\ 

= -Po{v±o - u±) + pi(v ±i - u±) - a— -— [R$mOa{y\\ -u\\)). (5) 

Here a is the internal column density of the shell, p is the volume density, v is the wind velocity, 
and u is the shell velocity. The independent variables are the spherical radius r, the spherical 
polar angle 9, and the time t, and R(9,t) is the radius of the shell. The quantities inside the 
bubble are labeled with the subscript "o" , while those outside the bubble bear the subscript "i," 
and quantities which are averaged over the width of the shell have no subscript. The velocity 
vectors have been decomposed into components perpendicular (e.g. u±) and parallel (e.g. uy) to 



-9- 



the shell. The quantity A, used for simplicity of notation, is defined by 



A = 



r 2 sin ( 



cos£ 



(6) 



R 



where £ is the angle between the radius vector and the local normal to the shell surface. The 
defining equation for £ is 



tan£ = — 



1 OR 
R~d9' 



(7) 



The first two terms on the right-hand side (rhs) of equation (6) represent the flux of mass through 
the inner and outer surfaces of the shell. The third term represents the change in surface density 
due to the change in surface area of the bubble, while the final term accounts for flow within the 
shell. 

The perpendicular and parallel components, respectively, of the momentum equation are 
given by 



a 



( dv± 
V dt 



l dt 



- [Po(v±o ~ U_l)(v± - V±)} 

+ [Pi(v±i ~ U±)(v_Li -V±)]+ Pi 



-{V\\ - 



"II 



a cos £ 
' R 



dv± 



-vii 1 + 



86 



(8) 



and 



V± di 



+ 



Po(v_L - U±)(V\\ - U|| 
pi(v ±i - U±)(v\\i -V\\) 



-{v\\ - 



"II 



(TCOSf 



dv\\ 



89 +V± { 1+ d6 



(9) 



where Pj is the thermal pressure causing bubble expansion, and we have assumed that P« S> P - 
In equations (8) and (9), the first two terms on the rhs represent the ram pressure balance. In 
equation (8) the third term represents the thermal pressure. The derivatives of and v± with 
respect to 6 in the final terms of equations (7) and (8) account for the ram pressure balance within 
the shell, and the term proportional to (1 + d^/d9) is due to the centrifugal force. 

The order of these equations may be reduced if one assumes that the flow is self-similar - i.e., 
that the expansion velocity is proportional to the spherical radius r. With this assumption, the 
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radius R may be written in terms of a self-similar radius L, using the expansion velocity at the 
pole, v p : R(6,t) = v p tL(9). For a RSG wind with constant mass flux M a and radial velocity v a , 
the density outside the bubble may be written 

Mgf(9) 

P ° = 4^#' (10) 

where f{9) is the asymmetry function, which is identically one for an isotropic wind. The 
self-similar surface density S(9) is defined by 

M a S{9) 



4:TtV a V p t 



(11) 



The parallel component of the gas velocity within the shell is given by 

v\\=v p U(9). (12) 
The parallel and perpendicular components of the shell velocity are defined by 

U|l = ~"^ Sm U± = ~dt C ° S ^ ' 

so that 

i£ 1 1 = — v p sm^L(6), u± = v p cos^L(9). (14) 
Because the BSG wind undergoes a strong shock, u± « v±. The velocities outside the bubble are 

v± = v a g{9) cos £, v\\ Q = -v a g(9) sin (15) 

where g{9) is the asymmetry of the RSG wind velocity. Because the RSG wind is assumed to have 
a constant mass flux, f{9) = l/g(9). We assume without loss of generality that /(0) = 1, and that 
L(0) = 1. Finally, the pressure inside the bubble is proportional to pv 2 . Given that the mass flux 
is constant, p ~ r -2 . Since v ~ r/t, we expect Pi ~ t~ 2 . Thus, we let P{ = F/t 2 , where F is a 
constant. 

With these substitutions, equations (5) and (7)-(9) become 

L' = -Ltanf, (16) 
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x ^-l) 2 -sece(A-l) 2 ] -1, (17) 

u> = f^-i)f^)-^d + 0, as) 



5 / V £/ + Lsin£ 

5' = S(tan£-cot#) _/(#)(£/ + L s i n £)- 2 

x f^-lV2C/ + (L + A„)sinO, (19) 



L 

where we have defined A = v a /v p , and = Xg(9). 

At the poles there is no flow within the shell, so that £(0) = U(0) = 0. In the limit — > 0, 
the self-similar form of the perpendicular component of the momentum equation (8) yields the 
equation 

* ' U (A-1) 2 . (20) 



A7rv a 



This has been used to eliminate F in equation (15). 

The value of v p is determined using the equation for bubble energy conservation. The rate of 
change of internal energy E is given by the work done by the fast wind, I, minus the work done to 
expand the bubble, dW/dt: 

dt dt y ' 

As is shown in DCB95, this yields the constraint on v p of 

1/9 /4«f\ 1/2 

[ J ^ L ) =°> ( 22 ) 

where Mf is the BSG mass loss rate, and x is the ratio of the volume of the bubble to the volume 
of a sphere with the same polar radius. 

The system of equations of solved as follows: A value of A is assumed. Equations (16)-(19) 
are then solved numerically to determine the self-similar functions L(6), S(9), U{6) and £(0). 
Given L(9) it is possible to calculate numerically x- Equation (22) then yields a new value of 
v p , and thus a new value of A. Iteration of this procedure may be used to find self-consistent 
energy-conserving models of the nebula surrounding SN87A. 
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The four variables v a , Vf, Mf and M a are parameters of the system which are varied to match 
observations of SN87A. The RSG wind density asymmetry function f{9) is determined using the 
wind-compressed disk model, described in §2.2, and the velocity asymmetry function g{9) is, as 
mentioned above, given by its inverse. 

Finally, we note that an implementation of the fifth-order Runge-Kutta algorithm was used 
to solve the WBB equations (16)- (19). A Newton's method was used to solve the additional 
constraint of equation (22). 

2.2. Wind-Compressed Disk Model 

The Wind-Compressed Disk Model (Bjorkman & Cassinelli 1993, hereafter BC) was developed 
to determine the structure of a rotating stellar wind. Their model applies whenever the external 
forces driving the wind are central {i.e., , radial) forces. BC found that, in the supersonic portion 
of the flow, the wind streamlines originating at high latitudes converge toward the equator. The 
reason for this behavior may be understood as a simple consequence of orbital dynamics. As the 
fluid leaves the surface of the star, it tends to orbit around the star until it is accelerated radially 
outward. If the initial outward acceleration is small compared to the rotation rate of the star, then 
the material orbits far enough around the star that it crosses the equator, where it meets material 
from the opposite hemisphere of the star. Since the flow velocities perpendicular to equator are 
supersonic, a pair of shocks form above and below the equator, and the shock compression of the 
material entering the equatorial region creates a dense disk. Thus we see that whether or not 
a disk forms depends on the rotation rate of the star, £1 = V ro t/V cr it = £1*/£Ik, where V ro t and 
V^ r it are the rotation and break-up speeds, and fi* and Qk are the stellar and Keplerian rotation 
rates. If the star rotates faster than a threshold value, > f2 t iu a dense, shock-compressed disk 
forms, with p ~ 100. Otherwise, there is only a mild density compression in the equator (typically 
p < 3-10). BC originally applied this model to investigate if the disks around Be stars could be 
produced by the WCD mechanism. The WCD model was later applied to several other classes of 
stars, including AGB stars, by Ignace, Cassinelli, & Bjorkman (1996, hereafter ICB). 
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To model the RSG wind, we follow ICB to obtain the flow velocity and density of the wind. 
ICB assume that the wind velocity of AGB stars can be described by the usual /3-law, 

v r (r, 9)=v + [ Voo (6 ) - v ] (1 - R*/rf , (23) 

where the initial velocity vq equals the sound speed, and is the wind terminal speed. The value 
of the velocity law exponent, /3, determines the acceleration of the wind. ICB assume that the 
winds of AGB stars are slowly accelerating, so they adopt (3 = 3 in their model. 

The wind density is determined by the continuity equation, which yields 

^- iW <24> 

where fi = cos 6, no = cosOq, and 9q is the initial co-latitude (at the stellar surface) of the 
streamline passing through the point (r, 9). Note that the factor (d/x/d/xo), given by equation 
(B3) of ICB, is the solid angle of the streamtube; this change in cross-sectional area produces 
the equatorial wind compression and is primarily responsible for the increased density near the 
equator. 

To evaluate the density, we must first determine 9q for the given location (r, 9) by solving 
numerically equation (B2) of ICB. However for our wind-blown bubble model, we are only 
interested in distances that are quite far away from the star. In the limit r — > oo, equation (A6) of 
ICB for the azimuthal deflection of the streamline becomes 

y rot sinfl ^ vq \W (I 1\ 

where B is the incomplete gamma function with the argument y m ax(#o) = 1 — vq/v 00 {9q). In the 
large-r limit, #o is given by the solution to the equation 

COS #o COS 0max(^o) = COS 9. (26) 

With this value for 9q, we then evaluate the wind compression factor (d^/d/io), given by equation 
(B3) of ICB, as well as the streamline mass loss rate and terminal speed, which are given by 

M(0o) = M a (1 - n sin 9 )t, (27) 
Voo(9 ) = v a (l -n sin # ) 7 , (28) 
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with 7 = 0.35 and f = -0.43 (see BC). 

Substituting equations (24), (27), and (|28|) into equation (0) and (0), we find the wind 



asymmetry functions 

g{9) = (l-n S m9 )\ (29) 
f(9) = (l-nsmeo)*^ (dti/d/io)- 1 . (30) 

These equations are valid for the wind. However, if the stellar rotation rate is above the disk 
formation threshold (Q > fith)i then we need the disk properties as well. Unfortunately, the WCD 
model per se cannot determine the disk density and velocities, so we make a simple estimate based 
on the numerical hydrodynamics simulations by Owocki, Cranmer & Blondin (1994). They found 
that the disk has an opening angle A9d ~ 3°, which agrees with observations of Be stars (Wood, 
Bjorkman & Bjorkman 1997), and that the disk terminal speed Vd ~ 0.2-0.3w a . However based on 
observations of the disks around Be and B[e] stars (Waters et al. 1988, Zickgraf et al. 1996), we 
adopt a slightly lower value Vd = 0.1v a . In addition to the geometry and flow speeds, Owocki et 
al. found that there is a stagnation point in the disk. Interior to the stagnation point, the disk 
material falls back onto the star, while exterior to the stagnation point the disk material flows 
outward. At low rotation rates, this recirculation is not dominant, so for simplicity we assume 
that all the wind material entering the disk flows outward. Using the continuity equation, the disk 
density may be estimated by 

p d = Mi (31) 
2irr 2 sin A9dVd 

The mass loss rate entering the disk, 

1 r^o ■ 

M d = - M(0 O ) d m , (32) 
^ Jo 

is determined by /j,q = cos 9q, the minimum co-latitude of all streamlines that enter the disk (see 
Fig. 10 of BC), which is given by the solution to 

C«) = ir/2. (33) 
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Using the disk terminal speed, O.lv, 



and equation (31), we find that the velocity and density 



asymmetry functions for the disk are 



9d 



0.1, 

/^(l-^sinflo)^ 
g d sin A6 d 



(34) 



fd 



(35) 



These equations are to be used only when O > O t h and \9 — ir/2\ < A6d; otherwise, one should 



As was discussed in §1, the accepted WBB model for SN87A treats the inner ring as the 
"waist" of the bipolar outflow and assumes that the outer rings lie somewhere on the surface of 
the lobes or represent lobe edges. In what follows we do not attempt to explain the outer rings, 
and primarily use the inner ring to constrain the parameters of our models. 

We use a number of observed features of SN87A for comparison with the models: The radius 
of the inner ring, w ~ 6.3 x 10 17 cm, sets the length-scale of our solutions. The density of the 
inner ring has been found to be at least ~ (1 — 2) x 10 3 cm -3 (Lundqvist & Sonneborn, 1997). 
Based on optical images, the width of the inner ring, 5r, is approximately one tenth of its radius. 
We use this to convert from column density to volume density: p m = a /Sr. The radial velocity of 
the inner ring is w ro ~ 8.3 km s _1 (Meaburn et al. 1995). The quantities v m and p m are the primary 
observations used to constrain the parameters of our models. We note that the polar expansion 
velocity may be given in terms of the velocity asymmetry function g{9) as v p = v m f{ir/2). In 
addition to t> ro and p m , we have calculated the ratio of the cylindrical radius of the lobes at their 
widest points to the radius of the inner ring (which we denote R' ). This is compared with the 
ratio of the cylindrical radius of the outer rings to that of the inner rings (which we denote R a ). 
The observed value of R Q is ~ 2. 



use equations (29) and (|30|), 



3. 



Results 



3.1. Fiducial Models 



As described in §2.2, the angular dependence of the RSG wind depends on the ratio 17 of 
the progenitor's rotation rate to the breakup rotation rate. The equator-to-pole density contrast 
p of the RSG wind is therefore proportional to Q. In what follows, we assume that the RSG 
has a radius of 740 R , a mass of 20 M , and an effective temperature of 3300 K. Given these 
values, the WCD model shows that when Q ~ 0.3, the equator-to-pole density contrast becomes 
significant, p ~ 25, even if a shock-bounded "disk" does not form. For Q ;> 0.3, a thin excretion 
disk forms with a very high density contrast of p ;> 700. Upon driving a BSG wind into these 
environments we find we can produce a bipolar wind-blown bubble consistent with the nebula 
around SN87A. In particular, we can recover the observed properties of the inner ring with two 
models, one with and one without a disk. The shape, expansion velocity and surface densities of 
these models are shown in Figs. 1, 2, and their density asymmetry functions f(9) are shown in 
Fig. 3. The initial parameters and the ring properties for the two models are given in Table 1. 
In both cases the expansion speed of the ring is equal to ~ 8.3 km s -1 . (We have also created a 
model with a ring expansion speed of ~ 10km s _1 , in accordance with the observations of Crotts 
et ai, 1995, and find that it requires, for the model with the WCD, an SI lower by 0.02; This does 
not significantly affect the results given by the Binary Merger Model of §4.) The ring density is 
within a factor of two of observed estimates (Table 1). 

It is important to note that we can produce successful models without resorting to anomalously 
low values for the speed of the RSG wind at the poles. As was discussed in §1, previous studies 
of SN87A have been forced to use velocities as low as v a ~ 5 km s _1 to recover the low ring 
expansion speed, v m . The hydrodynamics of Wind Compression naturally reduces v a at the 
equator, since mass conservation requires that the density is inversely proportional to the velocity. 
The reduced wind speeds at the equator allow us to use larger, more canonical values, of the RSG 
wind speed. 

The models with and without a disk differ most notably in the ratio R Q . Assuming that 
the outer rings lie at the widest point on the lobes, the observed value is R Q ~ 2. The model 
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with the disk has R Q = 4.9, whereas for the model without the disk, R Q = 2.6. The difference 
can be explained as follows: When a disk forms, its high density severely inhibits expansion of 
the ambient shock at the equator. A discontinuity in the shell radius L(6) occurs just above 
(below) the equator in these models. The BSG wind is able to expand rapidly along the the 
top and bottom edges of the disk producing a strongly "wasp-waisted," or pinched bubble. It is 
possible that such configurations are an artifact of our calculations. Numerical simulations have 
demonstrated that the lobes of energy conserving bubbles can experience considerable shaping due 
non-linear effects associated with the wind shock. If the wind shock is not spherical but assumes a 
prolate geometry, then post-shock BSG gas is focused towards the poles (Frank & Mellema, 1996). 
The higher ram pressures associated with this "shock focusing" will tend to elongate the bubble, 
reducing R Q (BL94). Frank & Mellema (1994) found that higher values of p produce stronger 
shock focusing. These effects will not be seen with a self-similar method which assumes the shell 
is driven by an isobaric "hot bubble". An isobaric bubble can not be produced with an aspherical 
wind shock. 

Figs. 1 and 2 and the ring properties given in Table 1 demonstrate the success of the unified 
GWBB and WCD models for SN87A. This represents one of the principal conclusions of our study. 
We take the cases represented in Figs. 1 and 2 as our fiducial models. 

3.2. Variation of Parameters 

We have varied the five input parameters to determine (1) the effect each has on the properties 
of the solutions, and (2) the degree to which the input parameters constrain the solutions. In 
particular, we are interested in the range of Q necessary to match observations, given that the 
RSG and BSG wind parameters are within the canonical range of values (e.g., v a ~ 20 km s -1 , 
400 £ Vf £ 600 km s~\ M a ~ 10~ 5 M yr _1 , M f ~ 10~ 7 M yr" 1 ). We find we can bracket the 
appropriate values of f2. Our results show that solutions outside the range 0.3 <J Q c <J 0.36 do 
not produce reasonable results for the bubble parameters. We therefore confine our discussions to 
cases within that range. 
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n 


V a 


M a v f 


Mf 


v w 




R Q 




km s _1 


1O" 5 M yr" 1 km s" 1 


1O~ 7 M yr" 1 


km s _1 


cm~ 3 




Fiducial Model 1 : Wind Compressed Zone 




0.3 


20 


2 400 


1 


8.30 


0.87 


2.6 


Fiducial Model 1 : Wind Compressed Disk 




0.36 


20 


1 600 


3 


8.36 


7.0 


4.9 


variation of parameters: 




0.33 


20 


2 400 


1 


7.02 


1.9 


3.2 


0.3 


25 


2 400 


1 


8.23 


0.73 


2.7 


0.3 


20 


1 400 


1 


10.5 


0.47 


2.5 


0.3 


20 


2 600 


1 


10.7 


0.94 


2.4 


0.3 


20 


2 400 


3 


11.7 


0.97 


2.4 



Table 1: Dependence of the properties of the models (inner ring expansion speed v w , inner ring 
number density n ro , and ratio of width of outer to inner rings, R a ) on the input parameters (fi 
ratio of RSG rotation rate to breakup rotation rate, RSG wind speed v a and mass loss rate M a , 
BSG wind speed Vf and mass loss rate Mf). 

As mentioned above, we find two models which provide agreement with the observed ring 
expansion speed: one with a disk, and one without. As Q increases, so does the pole-to-equator 
density contrast (and thus R a ). This results in higher equatorial densities, and thus smaller 
central ring expansion speeds. For a critical value of f2 = £l c such that 0.3 <^ J) c <^ 0.36, the wind 
asymmetry function develops a discontinuity with respect to 6, representing the presence of a disk. 
Models with a disk have considerably higher equatorial density, and must have correspondingly 
higher values of the BSG speed and mass loss rate if the same ring expansion speed is to be 
obtained (see Table 1). 

Increasing Vf or Mf results in a greater BSG wind momentum flux. This causes the central 
ring to expand at a higher rate, so v w is higher. This enables the bubble to sweep up more of the 
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RSG wind per unit time, resulting in a larger central ring density. 

Models with a higher value of v a possess a weaker RSG shock, resulting in lower equatorial 
expansion speeds and lower equatorial densities. When the RSG mass loss rate is decreased, v m 
increases, due to the smaller inertia of the RSG wind. Also, a smaller M a means less mass is 
swept-up by the expanding bubble; thus, the ring density is lower. 

Table 1 shows that the combined WBB/WCD model can be strongly constrained by input 
parameters and observations. 



4. Binary Merger Model 

As was seen in the previous section, for reasonable wind speeds and mass-loss rates, the RSG 
must have been rotating at least at ~ 0.3 of the breakup rotation rate. This is higher than can be 
expected for a single star (Fukuda, 1982; cf. Heger Sz Langer, 1998). Thus, our results imply the 
presence of a binary companion in order to spin- up the primary (Soker, 1998). In this section we 
present a simple model for a binary merger, and find the companion mass needed to achieve the 
rotation rates found in §3.2. Since we assume the companion to be of a much lower mass than 
the primary, it is expected that once the primary fills its Roche lobe, a dynamical mass transfer 
process will ensue (e.g., Rasio & Livio, 1996) with the subsequent spiraling-in of the companion, 
inside the common envelope. 

In order to obtain a simple lower limit estimate for the ratio of the mass of the companion to 
the mass of the primary, q = M p /M c , we assume that in the process of coalescence, the companion 
mass is added to the core of the primary, while all of its angular momentum serves to spin-up the 
primary's envelope (see Fig. 4). The initial angular momentum of the system is given by 



k 



G{M P M C ) 2 ^'' ' 



(36) 



M p + M c 

where a is the distance between the center of masses of the primary and companion. Given 
that the angular momentum after coalescence is primarily that of the envelope, the final angular 
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momentum is 

l f = £M e R 2 e n r , (37) 

where O r = QQk is the rotation rate (and Qk the breakup rotation rate), £ is the square of the 
gyration radius, R e is the radius, and M e = M e Q = M e i is the envelope mass, which is assumed to 
remain the same during coalescence. Setting Zj = If and simplifying, we find 

M^ 1 M C M P (M P + M c ) _1 / 2 (M e + M C \Y X I 2 = £ [~J (38) 

We assume that the final radius is approximately equal to the "volume equivalent" radius of the 
Roche lobe of the primary, Ry = aRy(q), where Ry may be approximated by the function 

~ _ 0.49<? 2 / 3 
V ~ 0.6g 2 / 3 + ln(l + gi/3) (39) 

(Eggleton, 1983), valid for all q. Using the assumption that M c <C M p , equation (38) may be 
further simplified, as 

(l + g)- 1 *^ 2 ^. (40) 

If we assume that the total mass is given by M c + M p = 2OM [so that 
q = M p /M c = (2OM - M c )/M c = 20/M c - 1] (Podsiadlowski, 1992), we may solve for 
the lower limit on the companion mass: 

M c ss 20M e £QR v (q). (41) 

For a giant star, £ ~ 0.1 (e.g., Pringle 1974). Equation ^ is a nonlinear equation for M c which 
may be solved numerically. For these values of [M c + M p ) and £ and for f2 ^ 0.5, Ry is close 
enough to constant that M c is approximately given by the linear function M c ~ 1.617. As a result, 
for = 0.3, the minimum companion mass is ~ 0.5 Mq. For the model with a disk, the minimum 
companion mass is predicted to be slightly higher, M c ~ 0.6 M . 

We must also consider the fact that material may also be ejected in the CE phase solely 
because of the the transfer of orbital energy into the envelope. As others have noted this is 
indeed a second possible mechanism for the development of the ring. The difficulty with using 
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such a scenario to link the shape of the nebula with the progenitor is that no simple way exits 
to connect the properties of the binary with the shape of the ejected slow wind and hence with 
the eventual shape of the bipolar nebula. Calculation of the ejection process requires expensive 
3-D calculations and these only provide information about relatively small scales (R < 10 14 cm, 
Sanquist et al. 1998). Thus while envelope ejection is alternative possibility which should be 
explored the advantage of the WCD model is that one can make definite statements about its 
consequences for the nebula via purely analytical methods. 

The dividing line between the WCD and envelope ejection models is determined by acE, the 
ratio of the envelope binding energy to the binary orbital energy which can be expressed as (de 
Kool 1990); 



acE = I^L = 4ffl / M ^ (42) 

AE orb ai M c M e - a f M p M c v ; 

where and aj are the initial and final orbital separations. If oqe > 1 then the entire envelope 
will not be lost in the CE phase. Consideration of this equation for parameters used in our models 
shows that if much of the orbital energy is deposited at radii > .01i? e then acE > 1 and we can 
expect the WCD to create a equatorial disk. It is possible that some mass is ejected in the early 
stages of the CE before the WCD begins. In fact, Soker (1998) recently proposed such a model to 
explain the presence of the two outer rings. 

We note also Podsiadlowski's (1992) suggestion that a binary merger is responsible for the 
transition from the RSG to BSC While our model is independent of the Podsiadlowski scenario it 
useful to compare the time scale for the two processes. Podsiadlowski quotes time scales for the 
merger process of a few 10 3 years. An upper limit on the time scale for the life of WCD is given 
simply by the dynamical time scale for the equatorial slow wind material, 



twcd < r dyn = *g!L » 10 4 (43) 

V a 5 km s 

The dense torus resulting from the WCD does not, however, have to be as large as the ring is 
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currently. It could have had a much smaller spatial extent (R ss 10 16 ) and still produce a nebula 
of the right size. Dwarkadas & Balick (1998), have shown that even a very small (R « 10 14 ) dense 
ring can hydro dynamically shape a bipolar nebula. Thus the time scales needed to create the BSG 
and the WCD appear close. 

We note also from Iben h Livio (1993) that "fcE, the ratio of the spin up time scale to the 
orbital decay time scale has the following proportionality 



1CE = ISSI. K P_, (44) 

Tdecay P 

where p is the local density the secondary star experiences as it spirals inward and p is the 
(considerably larger) average density for the RSG. Thus for an extended star like a RSG, jce < 1, 
and the envelope will be quickly spun up allowing the WCD mechanism to operate. 



5. Discussion and Conclusions 

We have developed a model for the bipolar nebula surrounding Supernova 1987A. As in 
previous investigations, we invoke the interaction between a fast isotropic Blue Supergiant Wind 
(BSG) and a slow aspherical (toroidal) Red Supergiant Wind. The novel aspect of our model is the 
the use of the Wind Compressed Disk (WCD) model to determine the geometrical properties of the 
RSG wind. Wind compression occurs when a star rotates fast enough to deflect wind streamlines 
towards the equator. The increased density and reduced velocity in the equatorial-zone RSG wind 
provides the constraint which shapes the final wind blown bubble (WBB). The degree of wind 
compression depends on the characteristics of the star. Thus our model allows the observable 
nebular properties to be directly linked to the unobservable properties of SN87A's progenitor. 

Our results show that the combined WCD/WBB model can recover the observed size, speed, 
density and gross morphology of SN87A's circumstellar nebula. Wind compression allows our 
models to recover the low expansion speeds of the ring without resorting to anomalously low 
values of the RSG wind speed. In addition, we have shown that our models are sensitive to initial 



-23- 



conditions allowing us to bracket the properties of the progenitor. In particular our models predict 
that the rotation rate of the progenitor must have been a significant fraction of the critical speed 
for break-up, Q ^ 0.3Qk- Since this is too large to be expected for a single RSG we infer that 
SN87A's progenitor was probably spun-up by a companion. Since no companion is visible now, 
we have developed a simple model for a binary merger. This model, along with the results of 
the WCD/WBB calculations, allows us to predict that SN87A had a companion with a mass of 
M c £ 0.5 M . 

Our value for M c clearly represents a lower limit. In addition to the fact that we have assumed 
an angular momentum deposition with 100% efficiency into the envelope, other effects also need 
to be considered. The shape of the bubble and speed of the ring are determined by the density 
contrast in the RSG wind. The formulation of the WCD model relies on an approximate method 
for tracing RSG wind streamlines. This ignores pressure effects. At large distances from the star 
it is likely that the wind compressed disk zone would experience some re-expansion which would 
weaken the density contrast. It is noteworthy that Garcia-Segura, Langer & MacLow (1997a) in 
their numerical WCD/WBB models for rj Carina found extreme stellar rotation rates ~ 0.9Qk 
were required to recover the correct morphology. The rapid evolution and high stellar temperature 
T* associated with r\ Car make it quite different from SN87A. However the Garcia-Segura models 
indicate that time-dependent hydrodynamic effects may change the properties of the disk. Thus, 
for SN87A, we expect larger values of would be required to produce high values of p and the 
correct inner-ring expansion speed. As equation (41) demonstrates, larger values of imply higher 
companion masses. 

In addition, we should note that the coalescence process is likely to produce some mass-loss 
which was not included in the calculation. A more accurate estimate of M c will therefore require 
both a complete common envelope evolution calculation as well as the use of numerical simulations 
of the bubble evolution. 

In spite of these uncertainties, our models provide a link between the properties of the 
circumstellar nebula and the progenitor star. Our results strongly indicate that SN1978A had a 
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binary companion which merged with the primary before the supernova explosion. A companion 
with a mass M c ~ 3 — 6M (as suggested by stellar evolution calculations; e.g., Podsiadlowski 
1992) is perfectly consistent with the dynamical considerations in the present work. 

We would like to thank Arlin Crotts, Larry Heifer and Peter Lundqvist for helpful discussions. 
Support for AF and TC was provided at the University of Rochester by NSF grant AST-9702484 
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Fig. 1. — The fiducial nebula model having no disk. The radius of the solid line shows the shape 
of the bubble, and is proportional to its expansion speed. The arrows indicate the magnitude 
and direction of the flow in the shell. The distance between the dashed line and the solid line is 
proportional to the surface density of the bubble. 
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Fig. 2. — The fiducial nebula model with a disk. As in Fig. 1, the radius of the solid line shows the 
shape of the bubble, and is proportional to its expansion speed. The arrows indicate the magnitude 
and direction of the flow in the shell. The distance between the dashed line and the solid line is 
proportional to the surface density of the bubble. The surface density jumps drastically near the 
equator, indicating the presence of a dense disk. 
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Fig. 3. — The density asymmetry functions for the RSG wind, for the models with (solid) and 
without (dashed) a wind-compressed disk. 



